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Suiniaary    In  the  recent  accounting  literature  there  has  been  much  in- 
terest in  the  construction  of  ARIMA  models  as  tools  for  forecasting 
corporate  earnings.   Our  examination  of  data  of  this  kind  indicates  that 
pov;er  transformations  of  the  observed  earnings  are  frequently  desirable. 
In  this  note  we  illustrate  the  model  building  procedure  for  two  quarterly 
series  of  earnings,  incorporating  into  the  usual  methodology  a  technique 
for  finding  an  initial  estimate  of  the  transformation  parameter,  v;hich 
in  turn  facilitates  the  identification  of  an  appropriate  ARIMA  structure. 

Key  words:   Time  series  model  building,   Transformation  of  data. 


1.   INTRODUCTION 
A  good  deal  of  recent  interest  in  the  Accounting  literature  has  fo- 
cussed  on  procedures  for  forecasting  corporation  earnings.   The  Financial 
Accounting  Standards  Board  (1977),  in  their  conceptual  framework  project, 
has  emphasized  the  importance  of  such  forecasts.    One  line  of  attack  of 
this  problem  has  been  through  the  construction  of  univariate  time  series 
models,  using  the  methodology  of  Box  and  Jenkins  (1970).   This  approach 
has  been  discussed  by,  for  example,  Foster  (1977),  Griffin  (1977)  and 
Lorek  (1979).   Much  of  the  research  has  concentrated  on  two  questions: 
do  corporate  earnings  streams  have  a  conmion  structure?   (that  is,  can 
one  find  a  single  model  from  the  general  autoregressive  integrated  moving 
average  class  which  predicts  v.'ell  for  a  wide  range  of  corporations?)  ;  and 
how  do  the  forecasts  from  time  series  models  compare  with  those  of  finan- 
cial analysts  and  management?   Some  discussion  on  the  latter  point  is 
contained  in  Abdel-khalik  and  Thompson  (1977-78),  Brown  and  Rozeff  (197S) 
and  Collins  and  Hopwood  (1980).   For  either  question,  it  is  clearly 
desirable  to  accumulate  evidence  of  the  time  series  behavior  of  a  large 
sample  of  corporate  earnings  series.   We  therefore  use  earnings  series 
as  examples  of  applications  of  the  technique  set  forth  in  the  present 
paper.   It  should  be  emphasized  that  our  primary  objective  is  not  to 
argue  for  or  present  evidence  favorable  to  individually  identified 
ARIKA  models  for  accounting  earnings,  but  rather  we  desire  to  demon- 
strate a  technique  applicable  to  a  v?ide  range  of  decision  oriented 
problems.   We  could  have  used  simulated  series,  but  we  wished  to  demon- 
strate the  technique  in  a  real  world  context. 
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This  paper  arose  from  the  study  of  a  large  collection  of  quarterly 
time  series  of  corporate  earnings  (Hopwood,  et.  al.,  1981).   It  quickly 
became  clear  that,  for  the  majority  of  series  in  our  sample,  there  was 
strong  evidence  of  the  desirability  of  a  data  transform.ation  to  induce 
homogeneity  of  error  variance,  although  this  point  had  not  previously 
been  noted  in  the  Accounting  literature,   0>ir  empirical  results  suggested 
that,  for  such  series,  a  transformation  will  generally  lead  to  forecasts 
of  improved  accuracy.   Accordingly,  we  considered  the  class  of  power 
transformations  of  Box  and  Cox  (196A).   Denoting  by  x   the  series  to  be 
analyzed,  the  class  of  models  considered  then  v.'as 

4,(B)$(B^)(1-B)'^(l-B^)°x^^^^  =  e(B)0(BSaj.  (1.1) 

where  a  denotes  white  noise  and  the  notation  is  that  of  Box  and  Jenkins 
(1970,  ch.  9),  with 

x^^^  =  (xj  -  1)/A  (A  /  0)  (1.2) 

Log  x^  (A  =  0) 

Along  these  lines  we  considered  four  possible  transformation 
strategies,  with  strategy  (d)  producing  the  most  accurate  forecasts: 

(a)  use  no  transformation 

(b)  use  the  logarithmic  transformation 

(c)  use  a  power  transformation  with  A  estimated  (jointly  with 
other  model  parameters)  by  the  maximum  likelihood  estimate  A. 

(d)  use  a  decision  rule  based  on  the  95  percent  confidence  in- 
terval for  A.   If  this  interval  contains  A  =  1  but  not  A  =  0, 
use  no  transformation.   If  the  interval  contains  A  =  0  but 


not  A  =  1,  use  the  logarithnic  transf oi-mation.   If  the  in- 
terval contains  neither  A  =  0  nor  A  =  1,  use  A  (-iointly 
estimated  with  other  model  parameters  by  maximuin  likelihood)  . 
If  the  interval  contains  both  A  =  0  and  A  =  1,  use  v;hichever 

is  closer  to  A. 

Under  the  usual  assumption  that  a   is  Gaussian  there  is  no  great 
difficulty  in  jointly  estimating  the  parameters  of  a  model  from  the  gen- 
eral class  (1.1).   However,  the  initial  specification  of  such  a  model  is 
com.plicated  by  the  fact  that  the  autocorrelation  structure  of  x  '   and 
its  differences  is  not  independent  of  the  choice  of  A.   Thus,  for  example, 
if  the  sample  autocorrelations  and  partial  autocorrelations  of  the  raw 
data  are  employed  in  the  usual  way  to  suggest  orders  for  the  autoregressive 
and  moving  average  operators  in  (1.1),  the  chosen  model  may  not  be  adequate 
to  describe  the  linear  properties  of  x    for  an  "appropriate"  A.   This 
point  is  established  theoretically  by  Granger  and  Newbold  (1976),  while 
a  numerical  example  in  Nelson  and  Granger  (1979)  shows  that  it  can  be 
practically  important.   In  our  study  we  found  that  the  problem  occurred 
in  about  10  percent  of  cases. 

In  the  next  section,  we  describe  an  elaboration  of  the  usual  model 
selection  procedure,  based  on  an  initial  estimate  of  the  transform.ation 
param.eter  A.  This  has  the  dual  advantages  of  providing  a  surer  initial 
identification  of  the  autoregressive-moving  average  structure  of  the  model 
for  X  and  of  yielding  the  preliminary  estimate  subsequently  required 
in  the  estimation  routine.  In  the  final  section  the  procedure  is  illus- 
trated through  the  modelling  of  two  corporate  earnings  series. 

2.   MODEL  IDENTIFICATION 
In  our  analyses  of  earnings  data  we  have  found  it  Vc-iluable  to 
base  model  selection  not  on  the  sample  autocorrelations  and  partial 
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autocorrelations  of  x  and  its  differences,  but  rather  on  those  of 

(X*) 
X     ,  where  \^   is  a  preliminary  estimate  of  A  in  (1.1)-(1.2).   This 

initial  estimate  is  obtained  by  fitting  autoregressions 


for  a  grid  of  values  A.   The  intention  is  to  obtain  a  reasonable  esti- 
mate of  the  variiince  of  a  ,  for  given  A,  in  (1.2)  through  the  usual  re- 

"2 
sidual  variance  a   ,  derived  from  fitting  (2.1).   To  achieve  this  we 
e ,  A 

choose  a  moderately  high  value  for  K,  the  order  of  the  approximating 

autoregression.   In  practice  setting  K  =  8  was  found  to  be  adequate  for 

2 
our  purposes   for  all  of  the  earnings  series  examined,  although  a  further 

elaboration  which  could  sometimes  be  useful  might  be  to  base  the  choice 

on  some  automatic  criterion  such  as  AIC  (Akaike  1974)  or  CAT  (Parzen 

1974).   The  initial  estimate  for  A  is  then  that  value  A*  whicli,  over 

the  chosen  grid,  maximizes 

g(A)  =  -l/2Kloga^   +  (A  -  1)  E   logx  (2.2) 

^'  t=l 


where  n  is  the  length  of  the  series  and  the  second  term  on  the  RHS  of 
(2.2)  is  the  logarithm  of  the  Jacobian  of  the  transformation  (1.2). 

Sample  autocorrelations  and  partial  autocorrelations  are  then  calcu- 

(X*) 
lated  for  x     and  its  appropriate  differences,  and  these  are  used  to 

identify  a  model.   This  modification  is  very  easily  incorporated  into 


existing  model  identification  routines  and,  since  (2.1)  can  be  ostirxited 
by  ordinary  least  squares,  the  additional  cost  in  terras  of  computing  time 
is  very  small. 

3.   MODELS  FOR  TOO  Ei"lRj\INC;K  SERIES 
The  data  analysed  in  this  section  consists  of  96  quarterly  earninp.s 
figures  for  two  corporations,  l^feyerhaeuser  Inc.  and  Freeport  Mineralr. . 
Fitting  autoregressions  (2.1)  with  K  fixed  at  8,  using  the  criterion  (2.2) 
and  searching  over  a  grid  of  width  0.05,  yielded  respective  initial  trans- 
formation parameter  estimates  a^'^  of  -0.19  and  -0.5A  for  the  two  serien. 
For  the  Ucyerhaeuser  data  the  sample  autocorrelations  of  the  transforn-ied 
series  indicated  that  a  single  non-seasonal  differencing  appeared  to  he 
sufficient  to  induce  stationarity.   Tlie  first  12  sample  autocorrelations 
and  partial  autocorrelations  of  the  differenced  series  are  shovTi  in  the 
upper  third  of  Table  1.   For  comparison,  the  middle  third  of  tlie  table 
gives  the  same  quantities  for  the  first  difference  of  the  untransformed 
series.   It  is  noticeable  that  the  magnitudes  and  patterns  of  the  two  sets 
of  sample  autocorrelations  are  quite  different,  particular] y  for  low  lags. 
From  the  values  for  A*  we  tentatively  identified  the  model 

(1  -  <J.E^)(1  -  E)x^^^  =  a  -  Qj  ^^)^^  (3.1) 

The  model  parameters  were  jointly  estimated  by  nuiT'erical  maximization  of 
the  exact  likelihood  function,  V7hich  is  readily  evaluated  by  incorporating 
in  an  obvious  way  the  algorithm  of  Ansley,  Splvey  and  TJrobleski  (1977)  into 
that  of  Ansley  (1979).   The  resulting  maximum  likelihood  estimates  (with 
estimated  standard  errors  in  brackets)  were 
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«>  =  0.38  (.10);    e  =  0.29  (.11);    >.  =  -0.28  (.18) 

The  bottom  third  of  table  1  shows  the  sample  riutocorrelations  and  partial 
autocorrelations  v/hen  the  data  is  transformed  according  to  A;  as  can  be 
seen,  these  are  virtually  identical  to  the  corresponding^,  quantities  for 
X* ,    indicating  that  the  initial  estimate  was  certainly  adequate  for  its 
required  purpose.   The  residual  autocorrelations  from  the  fitted  model  are 
shown  in  table  2.   In  general  these  are  not  sufficiently  high  to  cast 
serious  doubt  on  the  adequacy  of  the  fitted  nodal.   Given  the  figures  in 
table  1,  the  moderate  residual  autocorrelations  at  lags  7  and  9  are  not 
unexpected.   The  additional  elaboration  to  the  model  required  to  reduce 
their  magnitude  was  felt  to  be  unjustified. 

For  the  Freeport  Minerals  data  the  sample  autocorrelations  of  the 
initially  transformed  series  indicated  the  desirability  of  botli  a  non- 
seasonal  and  seasonal  differencing  factor  to  Induce  staticnarity .   Table  3 
gives  the  first  12  sample  autocorrelations  and  partial  autocorrelations  for 
the  appropriately  differenced  series.   The  most  important  difference  between 
these  and  the  corresponding  values  for  the  untrans formed  series  is  at  ]ag 
1,  where  in  the  latter  case  the  correlations  are  very  small.   The  auto- 
correlations for  the  transformed  data  suggest  a  model 

(1  -  B)(l  -  B^)x^^^  =   (1  -  eE)(l  -  OB^')a^  (3.2) 

The  maximum  likelihood  estimates  for  the  parameters  of  this  model  were 

e  =  0.23  (.11);    0  =  0.88  (.09);    A  =  -0.39  (.21) 

Again  it  can  be  seen  from  table  3  that  the  sample  autocorrelations  and 
partial  autocorrelations  when  A  is  used  differ  very  little  from  those  for 
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X*.  The  residual  autocorrelations  from  the  fitted  model  (3.2)  are  given 
in  table  4:  these  are  not  unduly  high,  suggestlnj;  that  the  cliosen  r.odel 
is  appropriate. 

A .   SUMMARY 
Our  experience  in  analyzing  time  series  of  earnings  data,  as  illus- 
trated in  the  previous  section, _^ indicates  that  very  often  a  power  trans- 
formation  is  desirable.   Certainly,  from  the  fitted  nodels,  the  hypothesis 
that  X  =  1  can  very  often  be  rejected  at  the  usual  significance  levels. 
VJe  also  found  that  the  modification  of  the  identification  procedure  dis- 
cussed in  section  2  is  typically  satisfactory  both  for  the  production  of 
sample  autocorrelations  and  partial  autocorrelations  on  which  a  reliable 
initial  model  identification  can  be  based  and  the  provision  of  an  initial 
estimate  of  the  transformation  parameter  for  use  in  a  maximum  likeJ.ihood 
estimation  routine. 
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FOOTNOTES 

For  example,  the  FASB  Statement  of  Financial  Accounting  Conceptis 
No.  1  (1978,  p.  ix)  stated,  "The  primary  focus  of  financial  reporting 
is  information  about  earnings  and  its  components."  The  reasoning  be- 
hind this  can  be  traced  to  the  tentative  conclusions  on  objectives  of 
Financial  Statements  of  Business  Enterprises  (FASB,  19  76)  where  the 
board  relied  on  four  propositions: 

(1)  ITie  primary  interest  of  the  investor  is  in  a  return  on  his 
investment  in  the  form  of  cash  flows  (p.  45). 

(2)  Earnings  measured  by  accrual  accounting  are  generally  thought 
to  be  the  most  relevant  indicator  of  an  enterprise's  cash- 
earning  ability  (p.  45). 

(3)  Fundamental  financial  analysis  focuses  on  earning  power  of 
an  enterprise  in  estimating  the  intrinsic  value  of  the 
stock  (p.  57) . 

(4)  Tlie  most  important  single  factor  determining  a  stock's  value 
is  now  held  to  be  the  indicated  average  future  earning 
power  (p.  57). 

Also,  forecasted  earnings  is  an  important  variable  in  research  involving 
cost  of  capital,  dividend  policy,  and  capital  markets  and  information 
content.   For  a  more  complete  discussion,  see  Foster  (1977,  p.  2). 

2 
Our  experimentation  with  different  values  of  k  revealed  that 

values  m.uch  smaller  than  8,  and  as  small  as  2,  produced  the  same  es- 
sential results. 
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Notes 

We  required  the  firms  to  be  listed  during  the  entire  sample 
period.  The  Center  for  Security  Price  Research  (CRSP)  monthly  tape 
was  used  to  select  NYSE  listed  firms.  A  firm  was  considered  listed 
if  it  had  monthly  stock  returns  available  for  the  entire  sample  period. 

2 
The  absolute  percentage  error  is  computed  as  the  average  of 

Actual  EPS  -  Predicted  EPS 


T  „„„         .  Since  this  error  metric  can  be  explosive 
Actual  EPS  '^ 

when  the  denominator  approaches  zero  we  truncated  errors  in  excess  of 
ten  to  a  value  of  ten.  This  operation  was  done  for  a  very  small  percent- 
age of  the  cases. 


